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Recently using Particle-In-Cell simulations i.e. in the kinetic plasma description Tsiklauri et al. 
and Genot et al. reported on a discovery of a new mechanism of parallel electric field generation, 
which results in electron acceleration. In this work we show that the parallel (to the uniform 
unperturbed magnetic field) electric field generation can be obtained in much simpler framework 
using ideal Magnetohydrodynamic (MHD) description, i.e. without resorting to complicated wave 
particle interaction effects such as ion polarisation drift and resulting space charge separation which 
seems to be an ultimate cause of the electron acceleration. In the ideal MHD the parallel (to the 
uniform unperturbed magnetic field) electric field appears due to fast magnetosonic waves which 
are generated by the interaction of weakly non-linear Alfven waves with the transverse density 
inhomogeneity. Further, in the context of the coronal heating problem a new two stage mechanism of 
the plasma heating is presented by putting emphasis, first, on the generation of parallel electric fields 
within ideal MHD description directly, rather than focusing on the enhanced dissipation mechanisms 
of the Alfven waves and, second, dissipation of these parallel electric fields via kinetic effects. It 
is shown that a single Alfven wave harmonic with frequency (u = 7 Hz), (which has longitudinal 
wavelength \a = 0.63 Mm for putative Alfven speed of 4328 km s~^) the generated parallel electric 
field could account for the 10% of the necessary coronal heating requirement. We conjecture that 
wide spectrum (lO""* — 10^ Hz) Alfven waves, based on observationally constrained spectrum, could 
provide necessary coronal heating requirement. It is also shown that the amplitude of generated 
parallel electric field exceeds the Dreicer electric field by about four orders of magnitude, which 
implies realisation of the run-away regime with the associated electron acceleration. 

PACS numbers: 96.60.-j;96.60.P-;96.60.pc;96.60.pf;96.50.Tf 



During the total solar eclipse of August 7, 1869, Hark- 
ness and Young discovered an emission line of feeble in- 
tensity in the green part of the spectrum of the corona. 
Similarly to the case of Helium discovered by Sir Nor- 
man Lockyer in 1868 it was mistakenly proposed that 
this line was due to an unknown element, provisionally 
named coronium. It was only in 1939 when Grotrian 
and Edlen correctly identified it as an emission produced 
by highly ionised Iron at temperatures of a few x 10^ 
K. Physical understanding of this high temperature in 
the solar corona is still a fundamental problem in astro- 
physics. None of the existing models can simultaneously 
account for all the observational and physical require- 
ments in order to explain what is known as the coronal 
heating problem. 

The second law of thermodynamics states that: It is 
not possible for heat to flow from a colder body to a 
warmer body without any work having been done to ac- 
complish this flow. Or energy will not flow spontaneously 
from a low temperature object to a higher temperature 
object. This seems to be in apparent contradiction with 
what we see in the solar atmosphere where temperature 
steeply rises from the photospheric boundary which is 
at T = 5785 K to a few x 10^ K i.e. an increase of a 
factor of 200. The transition region between the chro- 
mosphere (a thin layer above photosphere) and the hot 
corona over which the temperature increases dramati- 
cally is extremely thin (« 0.01 % of the Sun's diameter). 
This implies that some intense heat deposition occurs in 
the corona such that the second law of thermodynamics 
can be saved. 



The temperature structure of the solar corona is far 
from homogeneous Jll . The optically thin emission from 
the corona in soft X-rays or in the extreme ultra vio- 
let implies over-dense structures: so called coronal loops 
(closed magnetic structures) or plumes in polar regions 
(open magnetic structures) amongst others. In the 
corona the thermal pressure is generally much smaller 
than the magnetic pressure so that their ratio which is 
known as plasma-/? parameter is 1. This implies that 
in the corona substantial amount of energy is stored in 
the magnetic field. In turn, if one could devise a mecha- 
nism how this energy is released (dissipated) , then impor- 
tant clues to the solution of the coronal heating problem 
would be revealed. 

The coronal heating models are subdivided according 
to the possibilities of how the currents that are respon- 
sible for the plasma heating are dissipated: either by 
magnetic reconnection 0, S 3 ' Ohmic dissipation via 
current cascading 5, and viscous turbulence "T] in 
the case of so called DC models JU, or by Alfvenic res- 
onance, i.e. resonant absorption M [gl [iH [Tl[. phase 
mixing CIM III El E 111 13 MMllIlBl, and 
turbulence [23 in the case of AC models. An interesting 
alternative to all of the above models was explored in 
|2^. where coronal loops could be heated by the flow of 
solar wind plasma (plus other flows that may be present) 
across them by generating currents in a similar manner as 
a conventional Magnetohydrodynamic (MHD) generator. 

Historically all phase mixing studies have been per- 
formed in the Magnetohydrodynamic (MHD) approxima- 
tion, however, since the transverse scales in the Alfven 



2 



wave collapse progressively to zero, the MHD approxi- 
mation is inevitably violated. Thus, Refs.[^|23| studied 
the phase mixing effect in the kinetic regime, i.e. beyond 
a MHD approximation, where a new mechanism for the 
acceleration of electrons due to the generation of parallel 
electric field in the solar coronal context was discovered. 
This mechanism has important implications for various 
space and laboratory plasmas, e.g. the coronal heating 
problem and acceleration of the solar wind. It turns out 
that in the magnetospheric context similar parallel elec- 
tric field generation mechanism in the transversely inho- 
mogeneous plasmas was reported before [2^ |2^ (see also 
the comment letter and references therein). 

Contrary to the previous studies l^^, , here 

we use MHD description of the problem. Namely we 
solve numerically ideal, 2.5D, MHD equations in Carte- 
sian coordinates, with plasma /3 — 0.0001 starting from 
the following equilibrium configuration: A uniform mag- 
netic field Bq in the z— direction penetrates plasma with 
the density inhomogeneity across cc— direction, which is 
varying such that it increases from some reference back- 
ground value of po, which in our case was fixed at 
Pq — 2 X lO^pnip g cm^"^ (with molecular weight of 
p = 1.27 corresponding to the solar coronal conditions 
^Hi'^He^lGil JJ and irip being the proton mass), to bpa. 
Such density profile across the magnetic field has steep 
gradients with half-width of 3 Mm around x ~ ±10 Mm 
and is essentially flat elsewhere. Here "M" in units stands 
for mega i.e. 10^. Such a structure mimics e.g. foot- 
point of a large curvature radius solar coronal loop or 
a polar region plume with the ratio of density inhomo- 
geneity scale and the loop/plume radius of 0.3, which is 
a median value of the observed range 0.15 - 0.5 We 
use the following usual normalisation B^^y.z = BoB^ y^z, 
(x, y, z) = ao(a;, z), t = (ao/c^)t. Here Bq we fix at 100 
G, and hence dimensional Alfven speed, c\ = Bq/ ^/AiipQ 
turns out to be 4328 km s^^=:0.0144 c (c is the speed of 
light). The reference length oq was fixed at 1 Mm, i.e. 
dimensionless time of unity corresponds to 0.2311 s. We 
usually omit bar on top of normalised quantities, hence 
when simply numbers are quoted without units, in such 
circumstances we refer to dimensionless units as defined 
above. The dimensionless Alfven speed is normalised to 
c\. The simulation domain spans from — 40Mm to 40 
Mm in both x~ and z— directions with the density ramp 
as described above mimicking a footpoint fragment of a 
solar coronal loop or a polar region plume. Our initial 
equilibrium is depicted in Fig. 1. 

The initial conditions for the numerical simulation are 
By = Acoii{kz) and Vy = —CA{x)By at t = 0, which 
means that purely Alfvenic, linearly polarised, plane 
wave is launched travelling in the direction of positive 
z's. The rest of physical quantities, namely, Vx and Bx 
(which would be components of fast magnetosonic wave 
if the medium were totally homogeneous) and Vz and 
Bz (the analogs of slow magnetosonic wave) are initially 
set to zero. Plasma temperature is varied as inverse of 
density so that the total (thermal plus magnetic) pres- 



sure always remains constant. We fixed the amplitude of 
Alfven wave A at 0.05. This choice makes Alfven wave 
weakly non-linear. 

The models of coronal heating using wave dissipation 
were focusing on the mechanisms (e.g. phase mixing) 
that could enhance damping of the Alfven waves. How- 
ever for the coronal value of shear (Braginskii) viscosity, 
by which Alfven waves damp of about = 1 m^ s~^, 
typical dissipation lengths (e-fold decrease of Alfven wave 
amplitude over that lengths) are ~ 1000 Mm, only invok- 
ing somewhat ad hoc concept of so called enhanced resis- 
tivity can bring down the dissipation length to a required 
value of the order of hydrodynamic pressure scale height 
At ~ 50 Mm. In this light (observation of mostly un- 
damped Alfven waves 29] (see however Ref.'3^) and in- 
ability of classical (Braginskii) viscosity producing short 
enough dissipation length), it seems reasonable to focus 
rather on the generation of parallel electric fields which 
would guarantee plasma heating should the energy den- 
sity of the parallel electric fields be large enough. We 
put emphasis on parallel electric fields because in the di- 
rection parallel to the magnetic field electrons (and pro- 
tons) are not constrained. While in the perpendicular to 
the magnetic field direction particles are constrained be- 
cause of large classical conductivity cr = 6 x 10^^ s^^ (for 
T = 2MK corona) and inhibited momentum transport 
across the magnetic field. 

In the MHD limit there are no parallel electric fields 
associated with the Alfven wave. In the single fluid, ideal 
MHD (which is justified because of large cr) the parallel 
(to the uniform unperturbed magnetic field) electric field 
can be obtained from 



This means that in the considered system Ez can only 
be generated if the initial Alfven wave {Vy, By) is able to 
generate fast magnetosonic wave {Vx, Bx). 

In the recent past, in the context other than paral- 
lel electric field generation, Ref. investigated just 
such a possibility of growth of fast magnetosonic waves 
in a similar physical system. They used mostly analyti- 
cal approach and focused on the early stages of the sys- 
tem's evolution. Later, long term evolution of the fast 
magnetosonic wave generation was studied numerically 
in the case of harmonic [31 and Gaussian [16] Alfvenic 
initial perturbations. When fast magnetoacoustic pertur- 
bations are initially absent and the initial amplitude of 
the plane Alfven wave is small, the subsequent evolution 
of the wave, due to the difference in local Alfven speed 
across the x-coordinate, leads to the distortion of the 
wave front. This leads to the appearance of transverse 
(with respect to the applied magnetic field) gradients, 
which grow linearly with time. The main negative out- 
come of the studies with the harmonic 0| and Gaussian 
|16| Alfvenic initial perturbations was that the ampli- 
tudes of Vx and Bx after rapid initial growth actually 
tend to saturate due to the destructive wave interference 
effect [3. 
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We used fully non-linear MHD code Lare2d with 
the initial conditions as described above. As a self- 
consistency test we performed numerical simulations with 
moderate values of A: = 1 (in dimensional units this cor- 
responds to an Alfven wave with frequency {v — 0.7 Hz), 
i.e. longitudinal wave-numbers Xa = 6.3 Mm) and cor- 
roborated previous results of 0, E| . 

In Fig. (2) we show two snapshots of and (the 
latter was reconstructed using Eq.(l)) for the case of 
fc = 10 (in dimensional units this corresponds to an 
Alfven wave with frequency {v ^ 7 Hz), i.e. longitu- 
dinal wave- number = 0.63 Mm.). It can be seen 
that the fast magnetosonic wave (T4) and parallel elec- 
tric field {Ez) are both generated in the vicinity of the 
density gradients x ~ ±10, eventually filling the entire 
density ramp. This means the generated parallel electric 
fields are confined by the density gradients, i.e. the so- 
lar coronal loop (which the considered systems tries to 
mimic) after about 20 Alfven time scales becomes filled 
with the oscilla tory in time parallel electric fields. Also, 
as in '2l'| and '25*1, the generated parallel electric field 
is quite spiky, but more importantly it seems that large 
wave-numbers (and frequencies) i.e. short wavelength 
now are capable to significantly increase the amplitudes 
of generated both fast magnetosonic waves (Vx) and par- 
allel electric field E^ shown in Fig. (3). In Fig. (3) we 
plot the amplitudes of Vx = and E^ = E'^ which 
we define as the maxima of absolute values of the wave 
amplitudes along x ~ ±10 line (which essentially track 
the generated wave amplitudes in the strongest density 
gradient regions) at a given time step. This amplitude 
growth is beyond simple scaling (cf. [13 )■ In the 
considered case E^ now attains values of 0.001. To verify 
the convergence of the solution, we plot the results of the 
numerical run with the doubled (4000 x 4000) spatial res- 
olution. The match seems satisfactory, which validates 
the obtained results. One could argue that at first sight, 
the obtained spikes in both fast magnetosonic waves (Vx) 
and parallel electric field E^ do not contain much energy 
as they are localised strongly in the wavefront. However, 
one should realise that in Fig. (2) case of a single har- 
monic with 7 Hz is considered. If one considers a wide 
spectrum of Alfven waves (see discussion of this conjec- 
ture below), then the harmonics with different frequen- 
cies will contribute to the generation of what will then 
be a "forest" of such spikes (for each harmonic). Also, as 
the frequency of Alfven waves decreases than much more 
regular (no longer spiky) wave structures are observed for 
both fast magnetosonic waves (Vx) and parallel electric 
field Ez ■ To demonstrate this point we present intensity 
plots of Vx and E^ at 20 for the case of k = 1, v = 0.7 
Hz, Xa = 6.3 Mm in Fig. 4. In the latter case the am- 
plitude of Vx saturates at = 0.05^ = 0.0025 as one 
might expect from the weakly non-linear theory, while 
at 3 X 10"^ 

It has been know for decades [sj that the coronal en- 
ergy losses that need to be compensated by some addi- 
tional energy input, to keep solar corona to the observed 



temperatures, are as following (in units of erg cm~^ s"-'^): 
3 X 10^ for quiet Sun, 8 x 10^ for a coronal hole and lO"" 
for an active region. Ref.0] makes similar estimates for 
the heating flux per unit area (i.e. in erg cm~^ s~^): 

f„ = i.„A. = 5xlO»(^)'(^). (2) 

where Eh ~ 10~^ erg cm~^ s"^. This yields an estimate 
of Fh « 2 X 10^ erg cm~^ s^^ in an active region of the 
corona with a typical loop base electron number density 
of He = 2 X 10^ cm-3 and T = 1 MK. 

In order to make appropriate estimates we first note 
that the energy density associated with the parallel elec- 
tric field Ez is Ee = e£'^/(87r) erg cm~'^, where e is the 
dielectric permittivity of plasma. The latter can be cal- 
culated from e = {Awpc^)/ B'^ . This expression for e is 
different from the usual expression for the dielectric per- 
mittivity [s^: e — 1 + (47rpc^)/i3^ because the displace- 
ment current has been neglected in the above treatment, 
which is a usual assumption made in the MHD limit. 
At any rate for the coronal conditions (p = 2 x 10^/zmp 
g cm~'^, /i = 1.27, B = 100 Gauss) the second term 
(47rpc2)/B2 = 4.8048 x 10^ « e > 1. ^From Fig. (3) we 
gather that the parallel electric field amplitude attains 
value of « 0.001. In order to convert this to dimen- 
sional units we use — 4328 km s^^ and B = 100 G 
and Eq.(l) to obtain E^ « {c\B/c) x 0.001 = 0.0014 
statvolt cm~^ (in Gaussian units). Therefore the energy 
density associated with the parallel electric field Ez is 
Ee = (4.8048 x 10^) x 0.0014V(87r) = 3.7471 x 10"^ 
erg cm~'^. In order to get the heating fiux per unit area 
for a single harmonic with frequency 7 Hz, we multiply 
the latter expression by the Alfven speed of 4328 km s~^ 
(this is natural step because the fast magnetosonic waves 
{Vx and Bx) which propagate across the magnetic field 
and associated parallel electric fields {Ez) are generated 
on density gradients by the Alfven waves. Hence, the 
fiux is carried with the Alfven speed) to obtain 

Fe = Eec}]^ = 1.62 x 10^ [ergcm-^g-i] , (3) 

which is ~ 10 % of the coronal heating requirement esti- 
mate for the same parameters made above using Eq.(2). 
Note that the latter estimate is for a single harmonic with 
frequency 7 Hz (see discussion below for details when a 
wide spectrum of Alfven waves is considered). 

We now discuss how the energy stored in the generated 
parallel electric field is dissipated. For this purpose the 
parallel electric field behaviour at a given point in space 
(as a function of time) was studied. We found that Ez (in 
the point of strongest density gradient) is periodic (sign- 
changing) function which is a mixture of two harmonics 
with frequencies lo = CAk and u = 2c^fc. This can be 
qualitatively explained by the fact that Ez is calculated 
using Eq.(l) where Vy and By at fixed spatial point vary 
in time with frequencies lu = CAk, while the generated 
Vx and Bx do so with frequencies u = 2cAk 13]. Under 
influence of such periodic parallel electric fleld electrons 
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start oscillations and are accelerated (while ions perhaps 
are not due to their larger inertia). It should be empha- 
sised that because of the ideal MHD approximation used 
in this letter, strictly speaking the generated electric field 
cannot accelerate plasma particles or cause Ohmic heat- 
ing unless kinetic effects are invoked. Let us look at the 
Ohm's law for ideal MHD (Eq.(l)) in more detail denot- 
ing unperturbed by the waves physical quantities with 
subscript and the ones associated with the waves by 
a prime: Initial equilibrium implies i?o = Vq = with 
7^ 0. For the perturbed state (with Alfven waves 
{Vy and By) launched which generate fast magnetosonic 
waves {Vx and Bx) ) we have 

E' = -V' X (Bo + B')lc (4) 

Note that the projection of E' on full magnetic field (un- 
perturbed Bq plus the waves B') is zero by the definition 
of cross and scalar product: E' -{Bq+IS') /\{I3o^i3')\ = Q. 
Physically this means that in ideal MHD electric field 
cannot do any work as it is always perpendicular to the 
full (background plus wave) magnetic field. However the 
projection of i?' on unperturbed magnetic field Bq is 
clearly non-zero 

E'.^^E.^-f^.^ (5) 
\Bo\ c \B,\ 

which exactly coincides with Eq.(l) that was used to 
calculate E^ throughout this letter. Crucial next step 
needed to understand how the generated parallel (to the 
uniform unperturbed magnetic field) electric fields dissi- 
pate must invoke kinetic effects. In our two stage model, 
in the first stage bulk MHD motions (waves) are gener- 
ating the parallel electric fields, which as we saw cannot 
accelerate particles if we describe plasmas in the ideal 
MHD limit. Clearly we witnessed in works of Genot 
et al. and Tsiklauri et al. |23| that when identical 
system is modelled in kinetic regime particles are indeed 
accelerated with such parallel fields. Genot et al. 
claimed that electron acceleration is due to the polarisa- 
tion drift. In particular, they showed that once Alfven 
wave propagates on the transverse to the magnetic field 
density gradient, parallel electric field is generated due 
to charge separation caused by the polarisation drift as- 
sociated with the time varying electric field of the Alfven 
wave. Because polarisation drift speed is proportional 
to the mass of the particle, it is negligible for electrons, 
hence ions start to drift. This causes charge separation 
(the effect that is certainly beyond reach of MHD descrip- 
tion), which results in generation of the parallel electric 
fields, that in turn accelerates electrons. In the MHD 
consideration our parallel electric field is also time vary- 
ing, hence at the second (kinetic) stage the electron ac- 
celeration can proceed in the same manner through the 
ion polarisation drift and charge separation. The exact 
picture of the particle dynamics in this case can no longer 
be treated with MHD and kinetic description is more rel- 



evant. It should be mentioned that the frequencies con- 
sidered in kinetic regime v = 0.3a;ci/(27r) = 4.6 x 10^ 
Hz and this letter, 0.7 — 7 Hz, are different, but we clearly 
saw that the increase in frequency just yields enhance- 
ment of the parallel electric field generation. Various ef- 
fects of wave-particle interactions will rapidly damp the 
parallel electric fields on a time scale much shorter than 
MHD time scale. Genot et al. clearly demonstrated 
the role of nonlinearity and kinetic instabilities in the 
rapid conversion from initial low frequency electromag- 
netic regime, to a high frequency electrostatic one. They 
identified Buncman and weak beam plasma instabilities 
in their simulations (as they studied time evolution of 
the system for longer than Tsiklauri et al. |^ times). 
Fig. (11) from Genot et al. shows that wave energy 
is converted into particle energy on the times scales of 
lO'^w"/ ~ 4 Alfven periods. Perhaps no immediate com- 
parison is possible (because of the different frequencies 
involved), but still in our case Alfven period when sufh- 
cient energy is stored in the parallel electric field is 1/(7 
Hz) =0.14 s, i.e. wave energy is converted into the particle 
energy in quite short time. 

Our aim here was to show that parallel electric field 
generation is possible even within the MHD framework 
without resorting to plasma kinetics. Yet another impor- 
tant observation can be made by estimating the Dreicer 
electric field j32(|. Dreicer considered dynamics of elec- 
trons under action of two effects: the parallel electric 
field and friction between electrons and ions. He noted 
that the equation describing electron dynamics along the 
magnetic field can be written as 

rUeVd ^ eE ~ Vp^'-nieVd (6) 

where Vd is the electron drift velocity and Vp^^ is the elec- 
tron collision frequency and dot denotes time derivative. 
We note in passing that when Vd ^ vtuermai Eq.(6) in 
the steady state regime {d/dt — 0) allows to derive the 
expression for Spitzer resistivity. When Vd > vthermai 
then steady state solution may not apply. In this case 
when the right hand side of Eq.(6) is positive i.e. when 
eE > Vp^^meVd, we have electron acceleration. In fact 
this is so called run-away regime. In simple terms ac- 
celeration due to parallel electric field leads to an in- 
crease in Vd, in turn, this leads to decrease in Vp^^ because 
Vp^^ (X ^/v'^. In other words when electric field exceeds 
the critical value Ed = rieC'^ In A/(87r£gA:BT) (the Dre- 
icer electric field, which is quoted here in SI form), faster 
drift leads to decrease in electron-ion friction, which in 
turn results in even faster drift speed and hence run-away 
regime is reached. Putting in coronal values (rig — 2x 10^ 
cm~'^, T = 1 MK and In A = 17.75) in the expression for 
the Dreicer electric field wc obtain 0.0054 V m~^ which 
in Gaussian units is 1.8 x 10^^ statvolt cm~^. As can be 
seen from this estimate the values of parallel electric field 
obtained in this letter exceeds the Dreicer electric field 
by about four orders of magnitude! This guarantees that 
ran-away regime would take place leading to the electron 
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acceleration and fast conversion of the generated electric 
field energy in to heat. It should be noted, however, that 
in Eq.(6) E is constant, while our Ez in the point of 
strongest density gradient is time varying. Hence some 
modification from the Dreicer analysis is expected. 

It has been known for decades that the flux carried by 
Alfvcn waves in the solar corona is more than enough 
to heat it to the observed temperatures however 
linear MHD Alfven waves in homogeneous plasma do not 
possess parallel electric field component. This is one of 
the reasons why they are so difficult to dissipate. Hence 
the novelty of this study was to consider situation when 
parallel electric field generation is possible (weak non- 
linearity and transverse density inhomogeneity) . 

After comment paper by Mottez et al. [23| we came 
to realisation that electron acceleration seen in both se- 
ries of works 0, I25I is a non-resonant wave- 
particle interaction effect. In works by Tsiklauri et al. 
j2ll electron thermal speed was Vth,e = 0.1c while 
Alfven speed in the strongest density gradient regions was 
va = 0.16c, this unfortunate coincidence led us then to 
the conclusion that the electron acceleration by parallel 
electric fields was affected by the Landau resonance with 
the phase mixed Alfven wave. In works by [2fiL elec- 
tron thermal speed was Vth,e — 0.1c while Alfven speed 
was Va = 0.4c because they considered more strongly 
magnetised plasma applicable to Earth magnetospheric 
conditions. There were three main stages that lead to 
the formulation of the present model. 

(i) The realisation of the parallel electric field gener- 
ation (and particle acceleration) being a non-resonant 
wave-particle interaction effect, lead us to a question: 
could such parallel electric fields be generated in MHD 
approximation? 

(ii) Next we realised that indeed if one considers 
non-linear generation of fast magnetosonic waves in the 
plasma with transverse density inhomogeneity, then E = 
— (V X B)/c contains non-zero component parallel to the 
ambient magnetic field Ez = —iVxBy — VyBx)/c. 

(iii) From previous studies [la. we knew that the 
fast magnetosonic waves (14 and B^) did not grow to a 
substantial fraction of the Alfven wave amplitude. How- 
ever after reproducing old parameter regime (k=l, i.e fre- 
quency of 0.7 Hz), fortunately case of k=10, i.e frequency 
of 7 Hz was considered, which showed that fast magne- 
tosonic waves and in turn parallel electric field were more 
efficiently generated. 

We close this letter with the following four points: 

(i) In this work we showed that the parallel electric 
field generation re por ted in number of previous publica- 
tions 12^ which dealt with transversely in- 
homogeneous plasma can be explained in much simpler 
framework using MHD description, i.e. without resorting 
to complicated wave particle interaction effects. 

(ii) In the context of the coronal heating problem a 
new two stage mechanism of the plasma heating is pre- 
sented by putting emphasis, first, on the generation of 
parallel electric fields within ideal MHD description di- 



rectly, rather than focusing on the enhanced dissipation 
mechanisms of the Alfven waves and, second, dissipation 
of these parallel electric fields via kinetic effects. 

(iii) It is shown that a single Alfven wave harmonic 
with frequency = 7 Hz), (which has longitudinal wave- 
length Xa — 0.63 Mm for putative Alfven speed of 4328 
km s~^) the generated parallel electric field could account 
for the 10% of the necessary coronal heating requirement. 
We conjecture that the wide spectrum (10"^ — 10'^ Hz) 
Alfven waves, based on observationally constrained spec- 
trum, could provide necessary coronal heating require- 
ment. In this regard, it should mentioned that Alfven 
waves as observed in situ in the solar wind always appear 
to be propagating away from the Sun and it is therefore 
natural to assume a solar origin for these fiuctuations. 
However, the precise origin in the solar atmosphere of the 
hypothetical source spectrum for Alfvenic waves (turbu- 
lence) is unknown, given the impossibility of remote mag- 
netic field observations above the chromosphere-corona 
transition region [s^l ■ Studies of ion cyclotron resonance 
heating of the solar corona and high speed winds ex- 
ist which provide important spectroscopic constraints on 
the Alfven wave spectrum |35|. Although the spectrum 
can and is observed at distances of 0.3 AU, it can be 
then projected back at the base of corona using empiri- 
cal constraints, see e.g. top line in Fig. 5 from Cran- 
mer et al. |35|. Using the latter figure we can make 
the following estimates. Let us look at single harmonic, 
first. At frequency 7 Hz (used in our simulations) mag- 
netic energy of Alfvenic fluctuations is E^J^^'' « 10^ nT^ 
Hz"^ For single harmonic with = 7 Hz this gives for 
the energy density £;(7Hz) ^ „Ei'^^^^ /{8tt) 7 x lO'^ 
G^/(87r) « 2.8 X 10~* erg cm^'^. Note that surprisingly 
this semi-observational value is quite close to our theo- 
retical value of 3.7471 x 10"'* erg cm""^! As we saw above 
such single harmonic can provide approximately 10 % of 
the coronal heating requirement. Next let us look at how 
much energy density is stored in Alfven wave spectrum 
based on empirically guided top line in Fig. 5 from Cran- 
mer et al. [35^ . Their spectral energy density (which they 
call " power" ) is approximated by so called 1// spectrum, 
i.e. Ely = 0.6 X 10^ /v nT^ Hz~^. Which in proper en- 
ergy density units is E^, = 2.4 x 10~^/i^ erg cm~^ Hz 
Hence the flux carried by Alfven waves from say 10~^ Hz 
up to a frequency v would be 

F^"^ ^ r Ei,c\dy = 1.04 X 10^ ln(z//10-^) (7) 

JlO-'' Hz 

[ergcm~^s~^] . Based on the latter equation we deduce 
that the Alfven wave spectrum from v = 10~^ Hz up to 
about few xlO'^ Hz carries a flux that is nearly as much 
as the coronal heating requirement, Fh « 2 x 10^ erg 
cm~^ s~*, quoted above. We refrain from considering 
higher frequencies because for about 10"* Hz ions become 
resonant with circularly polarised Alfven waves and the 
dissipation proceeds through the Landau resonance - the 
well studied mechanism, but quite different from our non- 
resonant mechanism of parallel electric fleld generation. 
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There are several possibilities how this flux carried by 
Alfven waves (fluctuations), is dissipated. If we consider 
regime of frequencies up to 10^ Hz, ion cyclotron reso- 
nance condition is not met and hence dissipation would 
be dominated through the mechanism of parallel electric 
field dissipation formulated in this letter. However at 
this stage it is unclear how much energy could actually 
be dissipated. This is due to the fact that we only have 
two points 0.7 Hz and 7 Hz in our "theoretical spectrum" . 
As we saw a single Alfven wave harmonic with frequency 
7 Hz can dissipate enough heat to account for 10% of the 
coronal heating requirement. Therefore, we conjecture 
that wide spectrum (10^^ — 10'^ Hz) Alfven waves, based 
on observationally constrained spectrum, could provide 
necessary coronal heating requirement. More details on 



the issue of wide spectrum will be published elsewhere 




(iv) The obtained value of the generated parallel elec- 
tric field exceeds the Dreicer electric field by about four 
orders of magnitude, which implies realisation of the run- 
away regime with the associated electron acceleration. 
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FIG. 1: Dimensionless density {p{x)/po = 

[1 + 2 (tanh(a; + 10) + tanh(-a; + 10))]) (solid 
line) and Alfven speed, {ca{x) = 

l/v^l + 2 (tanh(.T + 10) + tanh(-x- + 10))) (dashed lino) 
profiles across the uniform unperturbed magnetic field 
(i.e. along .r-coordinate) which is used as an equilibrium 
configuration in our model of a footpoint of a solar coronal 
loop or a solar polar region plume. 
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FIG. 2: Snapshots of K and at t = 2 and 20 for the case of A: = 10, = 7 Hz, Xa = 0.63 Mm. 
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FIG. 3: Time evolution of the amplitudes of Vx = and 
Ez = E°,. Solid lines with stars represent solutions using the 
Lare2d code with 4000 x 4000 resolution, while dash-dotted 
lines with open symbols are the same but with 2000 x 2000 
resolution. Here k — 10, u — 7 Hz, Xa = 0.63 Mm. 



